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Abstract 



In this first paper we begin the application of variational methods to renor- 
malizable asymptotically free field theories, using the Gross-Neveu model as a 
laboratory. This variational method has been shown to lead to a numerically 
convergent sequence of approximations for the anharmonic oscillator. Here 
we perform a sample calculation in lowest orders, which show the superficially 
disastrous situation of variational calculations in quantum field theory, and 
how in the large- N limit all difficulties go away, as a warm up exercise for the 
finite-iV case. 
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I. INTRODUCTION. 



Many non-perturbative methods have been developed in recent years, trying to improve 
the asymptotic Taylor expansions (in terms of some coupling constant parameter) one usu- 
ally encounters in quantum field theories. A series of elegant papers in the seventies has 
explored and clarified various aspects of the large order behavior of perturbation theory 
in quantum mechanical and quantum field theoretical systems. In some cases this newly 
acquired knowledge was successfully used to obtain more accurate results from the asymp- 
totic perturbative series via Borel transforms. In other cases, especially in four-dimensional 
field theoretic models, the question of what the sum of a badly-behaved perturbation series 
means, remains unanswered |]J . 

Recently, there has been a considerable interest in variational-type expansions ( 0- [[UJ). 
Because of its conceptual simplicity and usually quite astonishing numerical accuracy, this 
type of expansion would be of particular importance in quantum field theoretical calcula- 
tions. Despite the distinctly intuitive nature of such a procedure, there has recently appeared 
strong evidence that optimized perturbation theory may indeed lead to a rigorously conver- 
gent series of approximants even in strong coupling cases. In particular, the convergence 
of this variational-like procedure has been rigorously established in the case of "zero" and 
"one" dimensional field theories [|lT| . 

In this paper, we begin the application to quantum field theory of the variational improve- 
ment of perturbation theory as described in [0. In the case of the anharmonic oscillator, 



this method gives a sequence of approximations which converges to correct answers as de- 



scribed in details in [13| , 14|] . Here, as a laboratory, using dimensional regularization, we 



study the 0(2N) Gross-Neveu model [15] in 2 + e dimensions. Exact results are known in 



2 dimensions for the mass gap |L6| and the vacuum energy density |L7| in terms of Ajjg for 
all N. 
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Before comparing variational calculations to these exact results, one must first under- 
stand how to handle infinities at e — > in a way compatible with perturbative renormal- 
ization, leaving only finite physical quantities expressed only in terms of finite physical 
renormalization group invariants of the effective variational approximation. This is the only 
sensible procedure if we want the physical characteristics of the theory, which have led to 
the specific choice of the effective variational approximation, to survive renormalization. 

In section [T| at e ^ 0, we compute the vacuum energy density and mass gap, according 
to the variationally improved perturbation theory of [P^-|14f in the few lowest orders, adding 



and subtracting a bare fermion mass, in closest analogy with the procedure which converges 
so remarkably in the anharmonic oscillator case. Not surprisingly, the procedure is in general 
highly singular as e — > 0, when the perturbative order is fixed. However, for N — > oo, the 
variational procedure gives the exact answers immediately when applied at first non trivial 
order. 

In section PTT| , using this remark and the fact that for iV — > oo the theory reduces to 
cactus diagrams which can be summed in closed (if implicit) form, we show how to remove 
the singularity of the procedure as s — > 0, obtaining finite, renormalization group invariant 
expressions for the vacuum energy and mass gap of the massless theory in term of the 
physical, finite, renormalization group invariant mass parameter of the effective massive 
theory, with respect to which we can then optimize. The optimum, which gives the exact 
known answers, turns out to be at zero, which is a physical consequence of the fact that one 
is actually solving the theory exactly. The procedure is sufficiently transparent to allow for 
a generalization to arbitrary N, which is done in the second paper. 
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II. LOW ORDER CALCULATIONS. 



In 2 + e dimensions, it is a straightforward exercise to compute the vacuum graphs of 
Fig. [I], which are the only ones which survives in the large N limit, using the Lagrangian 

N N 2 N _ 2 

c = %y, hWi + m o E + f- (E ^i) , (ii- 1 ) 

i=l i=l ^ i=l 

where m and go are bare parameters (in the following, we shall suppress the summation 
over the index i). We use the bare parameters for several reasons. First, from the remarkable 
results for the anharmonic oscillator we may conjecture that our procedure converges 

to the exact answer for fixed e < for any N. Second, the use of renormalized parameters 
introduces counterterms, which in effect mix perturbative orders and confuses matters as 
we are already perturbing in m®. It is thus preferable to perform as much of the calculation 
as possible with bare parameters, and only at the end relate them to the renormalized ones 
via the renormalization group through the standard procedure. The crucial requirement is 
of course that the final physical result remain finite for e —> 0. 

We give the value of the contribution of the graphs of Fig. [I] to the energy density E 
for all N : 



h-(,n) m ° +£jVr H) r 2 , r(-|) 

r 2 



+ (2N - l) 2 gt (1 + e) + 0(g«) 

= E^(m ) + g 2 Q E^(m ) + g* (m ) . (II.2) 



The straightforward application of the variational procedure at lowest order calls for the 
minimization with respect to m of the function 

(m ) = E® (m ) + gl E^ ( mo ) _ mQ ^^"° ) • (H.3) 
This gives the optimal value for the variational parameter 
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1 -( 2N - 1 )^$r9 2 o< {opt) , (0.4) 



and the corresponding optimum : 



rf- 



e s 



When iV goes to infinity, one can compute Eq exactly using the effective potential for the 
field a = In this large- iV limit, one finds that coincides with the exact Eq and mo 



with the exact mass gap for all values of e . This fact has already been noted in Ref. [E| and 



corresponds to the fact that for N = oo the theory is a set of iV free massive fermions, and 
the Oartree-Fock approximation is exact. We notice that for N = oo, g^N fixed, mo and 
E^ as given by Eqs. ( |II.4D and ( |0.5| ) are finite for e —>■ provided that the renormalized 



coupling g , given by 



% = - {N . 1)g T , ( n -6) 

7r e 



is kept fixed. Oowever, E^ and are separately infinite. 

The next order of the approximation calls for the minimization with respect to mo of 

E ( 2) (m ) = £ (0) (m ) + g 2 E^ (m ) + g' E^ (m ) 

dE<® (mo) 2 dE& ( mo ) 

- mo « fi-o m o a 

am am 

Just as for the anharmonic oscillator, and for the same reasons, Eq (mo) has an extremum 



at the same value of m [Eq. (|II.4|) 1 as Eq (m ), where it takes the same value [Eq. ( |0.5| )1, re- 



flecting the fact that in the large- A" limit, the variationally improved perturbative treatment 
gives the exact answer at each order. 

Although the behaviour of the large- A" limit is satisfactory, this sample calculation reveals 
two problems which must be solved before on can hope to get interesting results at finite A" 
or in more interesting theories. 



The first problem can be seen at lowest order, Eqs. ( [II. 4|) , ( [II. 5|) and ( [II. 6|) , which are 
true for all N: these equations do not give finite results for e — > at fixed finite N, due 
to the mismatch between the coefficient (N — 1) appearing in Eq. ( [11. 6| ) and the coefficient 
(iV — 1/2) appearing in Eq. ( [11. 4| ): the limits N oo and e — > do not commute at this 
stage, and it will be the essential content of this paper and the next one to make these two 
limits commute, i.e. to find how to reconcile in a general case a variational ansatz with the 
fine tuning of renormalization, materialized by Eq. ([II. 6|) and make the physics continuous 



as e crosses zero at fixed renormalized coupling constant g and mass scale //. 

The other problem appears at second order in Eq. ( pi.7|) , and is again related to the 
e — > limit. In this equation one finds that the quantity 

E<® (m ) - m + -m 2 v 2 u; II.8 

is finite as e —> although each term is separately divergent. Hence, in some sense, the 
lowest order in perturbation theory becomes irrelevant relative to the next two, each of 
order e^ 1 (taking into account the fact that g$ is of order e). This fact is independent of the 
large- N limit, and generalizes to higher orders of the calculation: at any finite order, for any 
N, the last two orders dominate all previous ones as e — > 0. This is quite unpleasant, as one 
would expect that all orders of perturbation theory contain at most a roughly equivalent 
amount of information. This is even more unpleasant in a perspective where one hopes to 
obtain accurate results from a low-order calculation! 

These two pathologies of the e — » limit will be handled in the next section in the 
large- N limit. 

Similar low-order calculations can be performed on the mass gap of the massive La- 
grangian ( [II. The contributing cactus diagrams up to order #q are on Fig. Q and, for all 
N, they give simply 

m F (m )=m + (2N - 1) g 2 m l +£ 



(47T) 1+ I 
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+ (2N - I) 2 g* mS+* (1 + s) + Ofog) 

= m? } + ^m? + g 4 Q mP . (11.9) 



Perturbing at order 1 in both and m , one has the mass gap 

drrip 
dm 



a (o) 

(1) / (0) 2 (1) OTflp o (1) 

my (m ) = my + g my - m — = g my . (11.10) 



This has no non trivial extremum in mo, but gives the exact large- N mass gap when the 
variational value of Eq. ( |11.4| ) is used for mo. Perturbing at order 2 in both g^ and mo gives 

(2) / \ 2 (i) 2 dmft 4 (2) m 
my (m ) = g my - 9 Q m -^- + g my . (H.ll) 



Again, this gives the exact large- iV mass gap when the value flll.4 ) for mo is used, but has 



a very singular behaviour if one tries instead to extremize it with respect to m . Indeed, in 
the limit e — > 0, the term g^m^ dominates, and it has no non trivial minimum. 

These pathologies extend to higher finite orders of the large- iV limit: plugging in the 
variational value ( |II.4|) for mo always gives the exact large-iV mass gap, but for e — ► 0, 
only the last order of perturbation theory survives, excluding any useful direct variational 
estimate of mp] just as for the vacuum energy, physics does not seem smooth in the e — > 
limit, and we shall also solve this problem in the next section. 



III. ALL ORDER CALCULATIONS AND RENORMALIZATION IN THE 

LARGE- N LIMIT. 



In this section, we fully exploit the fact that the large- N limit can be solved in closed 
(implicit) form and use further properties of the variationally improved perturbation theory 
discovered in one dimension (Ref. JHJ) to make renormalization transparent to our varia- 
tional method, and the physics continuous as e goes to zero. 
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For the vacuum energy density, the large- N limit amounts to retaining only the cactus 
vacuum diagrams. For the Lagrangian density of Eq. ( |11.1| ), this can be evaluated exactly 
and gives 

— ml N ml N 

= 2g g |1 - ^-^-f ry^r^ + 53 

+ ?T7TTW"* +E • (IIL12) 
2 + e (47t) + a 



where mp is given by 



m F = - ■ ( IIL13 ) 



1 - 



(2iV-l)ggm^r(-f) 
(4^) 1+ § 



Physically, in the usual treatment of the model, mp is the mass gap, and is determined 
from this last equation, which has a non trivial solution in the symmetric limit m — > 0. In 
the present variational method, we must take into account the perturbation in mo as shown 
in Eqs. ( |11.3| ) and ( |11.7| ). This is done rather easily in terms of contour integrals. Define a 
function f(x) by 

/(*0 = 1~ {2N T/lffi"^ x mg (1 - xf n ■ (IH.14) 

This function has power series expansions in both eft and x, beginning with 1. 

The n th order of perturbation theory for the vacuum energy density of the Lagrangian 
( gTIj ) is then: 

rnfnW \ 1 / dx _ nt . w Wig iV 

+ -L / ^_ /- , (1.1,5, 

2vk J x 2 + e (47r) + 2 

where the integration is counterclockwise on a small circle around the origin. 

It is of course straightforward to check that for n = 2 one precisely reproduces Eq. Qll.Tj ). 
It is also straightforward to check that 



v 0J = (111.16) 

orriQ 

at mo given by Eq. ( |II.4| ) and that the corresponding extremum of E^ n > is independent 
of n, establishing in a few lines that the variational result is independent of the order of 
perturbation theory in the large- N limit. Indeed, for 

l = (2N-l)^^g 2 ml , (111.17) 

one has simply 

f(x) = l-x . (III. 18) 



Furthermore in general one may rewrite the definition of /, Eq. (|IIL14|) as 



which shows that at fixed e < 0, / can be expanded in a power series of (1 — x). After a 
numerical exploration of the complex plane to make sure that no extra singularities lie in 
the way, one may distort the integration contour in Eq. ( |IH.15|) to run clockwise around the 
cut lying along the real positive axis and starting at x = 1. Actually, in Eq. ( |III.15|) , x = 1 
is an isolated pole, whose residue dominates its large n behaviour; this immediately shows 
that in the large- TV limit, at fixed e for any fixed m , E^ n \mo) goes to the exact answer 
Eq. (|II.5| ), in the limit of infinite order. This is reassuring: in the large N case, which is well 
under control, and where ordinary perturbation theory has a finite radius of convergence, 



the split of Ref. [12,13] between free and interacting terms is indeed independent of thq in 



the limit of infinite order, as intuition would suggest. 

However one could go one step further, and, as already noticed in the case of the an- 
harmonic oscillator in Ref. flT3"|, extract more structure from the limit of infinite order by 



rescaling mo with the order n. 

10 



After distortion of the contour it is clear that only the vicinity of x = 1 survives in the 
limit n — > oo, and one analyses this by the change of variable 

1 - x = - . (111.20) 

n 

Rescaling mo by introducing mo = m' n, m' kept fixed as n goes to infinity, one finds that 
[ n Eq (p;n,l5| ) has a limit E(m' ) given by 



E(m' ) = — I vdve v \- 
2m J 



/ 2 / 2 A7- 

///,, //;„ , A 



2g 2 f 1 (v) 2g 2i N-\ 

with /x defined simply by 

/ifr) - 1 ~ ^ 7/li!f ( ~ §) W /ifrP , (IIL22) 
(47rJ 2 

and the v integration contour lying counterclockwise around the negative real axis in the 
cut t>-plane. 

The function E(m' Q ) has the remarkable property, not shared by E^(m ), to allow for 
a smooth limit as e — > 0. This follows from the fact that Eq. (|III.22|) which defines f\ is 
continuous as e — > provided that g$ and mo follow their renormalization group behaviour. 
In this limit, we obtain simply 

T-,, /\ i i / t o \ m ' 2 r v dv e v N , _ 2 foil s /TTT 

E(m' )=E(m',g 2 ^) = — i —— — (1 + ^ 2 ) , (111.23) 



2z7T 7 /| 4vr v (AT - 1)^2 
with 

q 2 N m'v No 2 

f 2 (v) = l + ^ln— ^ - -^( lE + \nAn) 

7T /J,J 2 2l\ 

lE = 0.577215... , (111.24) 



where we have introduced the renormalized mass 
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m! 



™; = ^v7 • ( m - 25 ) 



i 

7T£ 



Eq. ( |111.23| ) involves only renormalized quantities, and is renormalization group invariant 
in the sense of the effective massive theory: 

[ V % + P{g) gj- g - i{g) m — } E{m', g\ /i) = , (111.26) 

with 

q 2 N q 2 N 

/% =-V~ ' -r(s) = — ■ ( m -27) 

Z7T 7T 

This follows trivially from the fact that E(m' , g 2 , fi) depends only on the bare parameters 
m' and g$. What is much less trivial is that we have arrived at a finite quantity for e — > 
(m ) in Eq. ( |111.15| ) is also renormalization group invariant but not finite in the limit 
e — > 0, while E(m' Q ) is both invariant and finite. 

One can go one step further by defining the dimensionless parameter 

m' it , 

m" = — (111.28) 

Ng 2 [ie 

and a function f 3 by 

/ 3 = ln(m%) - ln/ 3 . (111.29) 

m" is a finite pure number invariant by the renormalization group of the effective theory. 

Studying the function E(m") poses no problem. One finds that the function / 3 defined by 
Eq. ( |111.29| ) behaves as lnm" + 0(lnlnm") as m" goes to infinity, a typical renormalization 
groupimproved high energy behaviour. Furthermore, / 3 can be expanded in power series of 
v around v = 0, with radius of convergence 1/e, and with a cut extending from — 1/e ( a 
square root branch point ) to — oo. The behaviour of E(m") for large m" can be obtained 
by standard methods of contour integrals, and is of the form 
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m" 2 

E(m") ~ ^— . (111.30) 

This is the perturbative regime. For to" — > 0, E(m") converges exponentially to 

E(0) = -^fi 2 e'^ , (111.31) 

47T 

which is the only real extremum. Plotting the curve E(m") poses no problem with a work- 
station, and the result can be seen on Fig. [| Notice its remarkably smooth behaviour. 

The extremum at to" = is the exact value derived from the original calculation using 
the effective potential of the field a = iptp. This fact is consistent with the rescaling = 
m'o n as the order goes to infinity, and the fixed position of the extremum of E^ given by 
Eq. (fTJ). 



The same procedure can be applied to the calculation of the mass gap. From Eq. (pi.!3|) , 
the n th order of perturbation theory for the mass gap is 

(„) 1 £ dx _ n m 



m^> = — f> — x ~ n . (111.32) 

F 2m J x f K ' 

This reproduces the low order calculation ( |11.9| ); when mo is the value of Eq. ( [11. 4|) which 
extremises E^ n \ it is immediate that = m ; one may also go to infinite order, following 
the same rescaling of 1 — x and mo as for the vacuum energy density. In this infinite order 
limit, one obtains 

m F (m',g 2 ,fx) = — j —— , (111.33) 

a formula where the e — > limit has been taken like in the case of E(m', g 2 , fi), and contains 
the variational parameter m' in a finite renormalization group invariant way. 

mpfm', g 2 , fi) is actually a function of the renormalization invariant parameter to" only; 
TOi?(m") can be analyzed like E(m"), approaching exponentially the exact value at to" = 0, 
its only real extremum. It is plotted in Fig. f|. 
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IV. CONCLUSION. 



In this paper, with the Gross-Neveu model in the large N limit as example, we have 
studied some pathologies of the variational method applied to a renormalizable quantum 
field theory. On the two examples of the vacuum energy density and the mass gap, we have 
shown how a resummation to all orders can make these pathologies disappear, and make 
a variational method perfectly compatible with renormalization, i.e. one may replace the 
original theory ( the massless model ) by another one ( a massive model ), varying the 
parameters while keeping physics continuous as the space-time dimension is varied around 
its critical value. In the next paper, we extend the procedure to the finite-A case, where 
the mathematics of the renormalization group and the Feynman diagram stuctures are more 
generic. 
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Figure Captions. 



Figure 1. Vacuum energy graphs in the large- N limit. 
Figure 2. Mass gap graphs in the large- N limit. 

Figure 3. Variational vacuum energy E(m") / E(0) as a function of m" . 

Figure 4. Mass gap mp^m") /m^(0) as a function of the variational parameter m". 
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FIGURES 




FIG. 1. Vacuum energy graphs in the large- N limit. 
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FIG. 2. Mass gap graphs in the large- N limit. 
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FIG. 3. Variational vacuum energy E(m")/E(0) as a function of m" . 
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FIG. 4. Mass gap m j p(m")/m j p(0) as a function of the variational parameter m" . 
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